Macroscopic Quantum States and Quantum Phase Transition in Dicke Models of 

Arbitrary Atom-Number 
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The energy spectrum of Dicke Hamiltonians with and without the rotating wave approximation for 
arbitrary atom-number is obtained analytically with the variational method, in which the effective 
pseudo-spin Hamiltonian resulted from the expectation value in the boson-field coherent state is 
diagonalized by the spin-coherent-state transformation. In addition to the ground-state energy an 
excited macroscopic quantum-state is found corresponding to the south-and-north-pole gauges of 
the spin-coherent states respectively. Our results of ground-state energies in exact agreement with 
various approaches show that these models exhibit a zero-temperature quantum phase transition of 
second-order for any number of atoms, which however was commonly considered as a phenomenon 
of the thermodynamic limit with the atom-number tending to infinite. The critical behavior of 
geometric phase is analyzed, which displays no singularity at the critical point. 
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In this paper we analyze the energy spectrum of macro- 
scopic quantum-states (MQS) in relation with the quan- 
tum phase transition (QPT) for N two-level atoms inter- 
acting with a single bosonic mode. The system is known 
as the Dicke model (DM) [l], while the MQS refers to bo- 
son and spin coherent-states. Although the model itself 
is quite simple, it displays a rich variety of the unique as- 
pects of quantum theory and has become a paradigmatic 
example of collective quantum behaviors. Thus DM is 
of fundamental importance in quantum optics and any 
systems involving the interaction between matter and 
light, for example, the superconducting circuits coupled 
to a resonator. Since the atoms are confined in a small 
container compared to the wavelength of radiation-field 
the model is reduced to a quantum mechanical problem, 
which can be solved exactly. At zero temperature, the 
system of TV interacting atoms via the bosonic field can 
exhibit a transition to a superradiant phase called the 
QPT, which describes a structural change in the proper- 
ties of the ground-state energy spectrum [2j at the critical 
value of coupling parameter, where both the numbers of 
field-bosons and atoms in excited states exhibit an abrupt 
increase from zero. The QPT has become a fundamental 
way to explore dynamic properties of quantum correla- 
tions in many-body physics and thus is of central interest 
in the realm of quantum information theory ^3^ as well. 
The DM, which exhibits a second-order phase transition 
from the normalphase to the superradiant phase pre- 
dicted long ago [j, Q , has been studied extensively be- 
cause of its broad application range such as the ground- 
state entanglement and corresponding finite-size behav- 
ior. The critical exponents and finite-size corrections are 
obtained [6] by a Holstein-Primakoff (HP) approach Q, 
which is the most popular method [s^ in the investigation 
of QPT of DM. 

It is believed that the QPT can take place only if the 
collective atom-photon coupling strength is the same or- 
der of the atom level-space, which was considered as a 
challenging transition-condition. In the strongly coupled 
regime of cavity quantum electrodynamics this condition. 



however, is shown to be accessible by controlling the 
pump laser power Q. For a Bose-Einstein condensate 
in a high-finesse optical cavity, the energy space of two 
levels can be adjusted to be small enough and the QPT, 
namely the superradiance transition, has been observed 
experimentally [lol - [l^ . 

The HP representation that converts the Hamiltonian 
into a two-mode bosonic problem is the starting point for 
the most theoretical analysis of the QPT in the thermo- 
dynamic limit {N — > oo), which leads to solvable Hamil- 
tonians by neglecting terms from expansions of the HP 
square roots [8]. In this paper we present an alternative 
procedure to study the energy spectrum of MQSs with 
the standard variational method. An effective pseudo- 
spin Hamiltonian in the energy functional is diagonal- 
ized in terms of the spin coherent state (SCS) transfor- 
mation [l3|, which is the key point to lead to the 
new observation. Both the energy-spectra of two MQSs 
obtained from the south-and-north-pole gauges of SCSs 
are found to have the behavior of QPT, where only one 
is the ground state. In the MQSs the QPT regarded 
as a drastic change of ground-state energy spectrum ex- 
ists for any number of atoms since the SCS is valid for 
any spin- value including the spin-i corresponding to the 
one-atom case and the thermodynamic-limit approxima- 
tion resulted from the HP-transformation is avoided. The 
critical behaviors of geometric phase (GP) are also ana- 
lyzed. 

Dicke Hamiltonian and spin- coherent- state represen- 
tation Dicke Hamiltonian (DH) describing the interac- 
tion of TV two-level atoms of level-space VI with a single 
bosonic mode of frequency uj, can be generally written as 
H = Hq + Hsb, where Hq — uja^a -\- ftJz,B,nd the spin- 
boson coupling part is defined as 



Hsb — 



J'x {a) 



(1) 



The operators a and a) denote the one-mode annihila- 
tion and creation boson (photon)-operators respectively, 
Jz is the atomic relative-population operator, and J± 
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= (JxiiJy) are the atomic transition operators. Various 
approximations have been developed, which are suitable 
for particular ranges of parameters. The rotating-wave 
approximation (RWA), which has broad application in 
quantum optics, is based on the assumption of near res- 
onance and relatively weak coupling between the boson- 
field and atoms [T5| . The spin-boson coupling part be- 
comes 



h: 



2VN 



[Jx{a + a^) + iJy{a — a^)] 



(2) 



in the RWA. 

We are going to determine the energy-values of the 
MQSs as a function of the coupling constant g. To this 
end we use a variational method with the trail state being 
a direct product of SCS and Weyl (boson coherent) state 
[Hi IV'i) = ± n) where the boson coherent-state is 
defined by a\a) = a\a) and the SU{2) SCS is 



J • n| ± n) = ±s| ± n) 



(3) 



with s — Y being the total pseudo-spin value. The unit 
vector 



n =(sin 9 cos cj), sin 9 sin (j>, cos 9) 



(4) 



expands a Bloch sphere. The states | ± n) defined 
as an eigenstate of operator J • n with eigenvalues ±s 
are called the SCSs of north-and-south-pole gauges re- 
spectively, in which the minimum uncertainty ^KJx)] — 
((AJ^)2)5((AJy)2)5 is satisfied, where (J,) = (±n| J^|± 
n). It has been shown that the expectation value of 
spin in the SCS satisfies the exactly same Bloch equa- 
tion as that of a classical magnetization- vector in exter- 
nal fields 01, so that we refer it as a MQS. The SCSs 
can be generated from the extreme states of Jz, that 
Jz\s,±s) = ±s\s,±s), by a rotation | ± n) = i?|s,±s) 
with the unitary operator obviously being 



(5) 



where the unit- vector m in the x-y plane is perpendicular 
to the plane expanded by the z-axis and unit vector n. 
The energy expectation value in the trail state reads 



E±{a) = (^±|i?|^±) = {±n\Hes{a)\ ± n) 



(6) 



The effective spin-Hamiltonian H^s (a) resulted from the 
bosonic mean-field treatment contains the spin operators 
only, however, with the complex bosonic-field parame- 
ter a to be determined by the variation-principle which 
minimizes the energy functional E± (a) . The key point 
of our method is that we can properly choose the unitary 
operator R defined by Eq.® such that the SCSs | ± n) 
become eigenstates of the effective spin-Hamiltonian, i.e., 
Hes (a) I ± n) = E± (a) | ± n) , which describes N atoms 
in an effective, classical light-field parameterized by a. 
In this stage the energy functional is obtained exactly 
compared with the HP approach, where the square root 



is expanded approximately with large atom-number N. 
The ground state energy can be obtained by the standard 
variational method with respect to the complex param- 



eter a, 



da 



or 



OE+ja) 



0, in which only one 



equation can be satisfied giving rise to the ground-state 
energy-spectrum. For the general unit-vector n{6, (/>) in 
Eq.Q, the vector m to be determined is obviously 



:(sin0, — coS(/>, 0). 



(7) 



It has been demonstrated by the HP representation 
that in the thermodynamic limit (N — ^ oo) the sys- 
tem undergoes a QPT at a critical coupling gc, at which 
point the system changes from a largely unexcited normal 
phase to a superradiant one, where all atoms absorb and 
emit bosons collectively and coherently This QPT 
has recently been observed by employing a superfluid gas 
in an optical cavity pX)| . The theoretical description is 
based on HP realization of the quasi-spin operators by the 
infinite expansion-series of additional bosons. The infi- 
nite series are truncated up to the first or second order 
with the large N assumption, i.e. the number of excited 
atoms is small compared with the total number of atoms 
B B [3 ■ It is also shown that the number of photons is 
divergent [1, [l^ at the critical value of coupling constant 
gc- We in the present paper provide the analytical results 
of the energy spectrum of MQSs, the occupation-number 
of atom levels and the expectation value of photon num- 
bers in terms of the SCS transformation. The critical 
points of QPT are found explicitly. 

Energy spectrum of macroscopic quantum states and 
quantum phase transition For the DM with spin-boson 
coupling Eq.(IT]), the effective spin Hamiltonian is 



(8) 



where a = u + iv, and Hs{a) = VtJz + '^Jx denotes the 
effective spin Hamiltonian with the constant excluded, 
which can be written as Hs{a) = r{cos9Jz + sm9Jx) 

with the parameters determined as r = y^fJ^ -I- (^-^u)^, 

cos 9 — —, and sin 9 — ^g"/^^ jn this case the rotation 
operator Eq.(l5]) has a simple form with the unit vector 
m = (0,-1,0) seen from Eq.©. The energy eigenvalue 
of Hs{a) is Es{a) = sr and the total energy functional is 



E^{a) = u;{u' + v^)±^Jn^ + i^n)^. (9) 



The ground-state energy is the minimum of energy func- 
tional (a) for the SCS of south-pole gauge only. From 
— 0, ^"^a = we have v — 0, and m = or 



Nn^ fg 



1 



(10) 



which is photon-number in the superradiant phase. The 
parameter 



9c 



(11) 
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is just the critical point of phase transition and is the 
same value obtained with the HB approach [8]. Then 
the ground-state energy is 

f-f^, (9^ 9c) , , 

which is also obtained in Refs. [1, [l^, [l^, however is 
just an approximation of thermodynamic limit in the HP- 
transformation approach. We demonstrate that the ap- 
proximation resulted from the HP transformation itself 
can be avoided with the SCS representation. Below gc the 
system is in normal phase but the superradiation phase 
above the critical point, while the QPT emerges for any 
number of atoms in our formulation and it is no need 
at all for the thermodynamic limit. The smooth QPT 
is of the second order seen from the continuation of the 
first-order derivative -^\g=g^_^, = -^|g=g^„ = at the 
critical point (see Fig. [T]). The one half of occupation- 
number difference between two levels of the atoms can 
be evaluated by 

{J.} = {-s\l\-s) = \~f^^^ ^f^-^'l (13) 
[ ^Tp-, (9 > 9c) , 

where = B) J^R = cosdJz — sm9Jx- It is interesting 
to see the energy spectrum of excited MQS resulted from 
the SCS of the north-pole gauge that 

(9^ 9c) 

^-^l^(^-i^), i9>9c), ^''^ 

which although cannot be occupied at zero temperature 
exhibits the behavior of QPT at the same critical point 
gc- This excited MQS may be of importance in laser 
physics and technology. The energy spectra and the 

occupation-probability difference as a function of 
coupling constant g are plotted in Fig. [1] For the one- 
atom case s = i, (Jz) is nothing but the one half of 
occupation-probability difference between the two levels. 

Rotating wave approximation In quantum optics the 
DH is usually studied in the RWA for small value of cou- 
pling, in which the counter-rotating terms J+ and aJ_ 
are neglected so that the DH becomes integrable. The ef- 
fective spin Hamiltonian for the RWA is seen from Eq.([2]) 
as 

= nj, + -i= {J^U - JyV) , (15) 

V N 

which can be rewritten as the standard J • n form, 
i.e., i/J = r [cos^Jz -I- sin0 (cos^Ja; -I- sin^Jy)], where 
COS0 = -A, sine!) = -p^, cos9 — — , aind = M=-, and 

r = y^fl'^ + g'^{u'^ + v"^) /N . The unitary operator with 
the vector m given in Eq.© can transfer the state |s, ±s) 
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FIG. 1: (Color online) The coupling-constant dependence 
of energy-spectra of MQSs E^/N, E+/N (dashed line) and 
the occupation-probability difference between two atom-levels 
(Jz) /N (inset) with ui = Q. = 1 and the critical point Qc — 1. 



into the eigenstate of Hg, and the energy functional in 
this case becomes 

EUa)^.iu^ + v^)±^f^7^!^^. (16) 
From ^^Q^"^ = 0, and '^^g^"^ = we have u = v — oi 

where the critical point is seen to be 

gl - 2V^, (18) 

which is two times of the value gc in the DH because 
of neglecting the counter-rotating terms and is in agree- 
ment with the result obtained in Ref. 0, where it was 
found that the phase transition does not exist below this 
value. The same critical point is also obtained in Ref. 
[20| with HP transformation. The ground-state energy 
spectrum E'[_ has the exactly same form as Eg. p^ for 
the DH, however, with the critical point replaced by 
and so does the spectrum of excited MQS E!^. The 
occupation-number difference (Jz) is also formerly the 
same as Ea. ([T^ . 

Critical behavior of geometric phase Considerable un- 
derstanding of the time-evolution of quantum states has 
been achieved after Berry's discovery of its geometric fea- 
ture known as the Berry Phase. An apparently unrelated 
area in connection with CP, which has attracted consid- 
erable attention recently, is the study of QPT in many- 
body systems 20-22]. QPT is accompanied by a drastic 
change in the ground state at critical points, which is 
associated with a non-analytical behavior of the ground- 
state and is expected to be reflected in the CP of the wave 
function. In the DM with various modifications the CP 
is generated by a rotation with photon-number operator 
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FIG. 2: (Color online) The coupling-constant dependence of 
GP {"f/N) and its derivative d'yj{Ndg) (inset) with u = Q, = 1 
and the critical point = 1- 



as in the quantum optics such that [20|, [2 

7--» r {E^\U\^)^U{4>)\E^)d(p 
Jo ocj) 

f 0, ^ (5 < .9c 



27r|a| 



23^ 



(19) 



where the unitary transformation-operator is seen to be 
U{(j)) = e^"^""^. The hnear TV-dependence of GP (or 
photon-number) remains in the superradiant phase hav- 
ing no particular divergence at the critical point . The 
first-order derivative of the GP is obviously 




(3 ^ 9c) 
{9 > 9c) ■ 



(20) 



Besides the factor-2 difference of critical point {g^ = 2gc) 
the photon number in the RWA case given by Eq. pT)) 
has an extra factor-4 compared with the photon-number 
formula Eq.ldU]) for the DH. Fig. [5] shows the GP (7/A^) 
and its derivative {^j^) (inset) as functions of g for the 
DH. Unlike the energy spectrum, the GP exhibits a sharp 
change at the critical point with discontinuous first-order 



derivative. The scaling behavior of the GP at the critical 
point can be found as 



^1 



\9-9cl 



(21) 



which agrees with the result of HP representation in the 
thermodynamic limit [20] and shows no singularity at the 
critical point confirming the conclusion in Ref. [l6|. 

Conclusion The MQSs are useful in the description of 
QPT, which is, as a matter of fact, a macroscopic quan- 
tum phenomenon. The effective pseudo-spin Hamilto- 
nian can be diagonalized by the SCS transformation and 
the energy spectra of ground-state along with an excited 
MQS for any number of atoms are obtained, which pos- 
sess the exactly same form of formulas for both DHs with 
and without the RWA expressed in terms of the critical- 
point parameter, while the critical points themselves have 
a factor-2 difference. The SGS technique universal for 
various Hamiltonians of Dicke-type, is valid for strong 
spin-boson coupling while has no restriction on the num- 
ber of atoms. Thus we conclude that the abrupt change of 
ground-state energy spectrum regarded as QPT is an ef- 
fect of strong light-field, which is independent of the atom 
number in DHs. At the critical point of phase transition, 
the expectation values of photon and atom numbers in 
excited states exhibit an drastic increase from zero pro- 
portional to the total number of atoms N and remain the 
same in the superradiant phase. There is no singular be- 
havior at the critical point as demonstrated in Ref. (l6j . 
where the tensorial product of spin-coherent and Weyl 
states firstly used to analyze the ground-state energy of 
DM is shown to lead to an energy surface of the DH with 
four variables and two parameters. We adopt the SCS 
transformation to obtain the analytic energy spectra of 
ground state and an excited MQS as well for both the 
Hamiltonians with and without the RWA. 
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